Charge density wave and superconducting dome in TiSe2 from electron-phonon 

interaction. 
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At low temperature TiSe2 undergoes a charge density wave instability. Superconductivity is sta- 
bilized either by pressure or by Cu intercalation. We show that the pressure phase diagram of TiSe2 
is well described by first-principles calculations. At pressures smaller than 4 GPa charge density 
wave ordering occurs, in agreement with experiments. At larger pressures the disappearing of the 
charge density wave is due to a stiffening of the short-range force-constants and not to the variation 
of nesting with pressure. Finally we show that the behavior of Tc as a function of pressure is entirely 
determined by the electron-phonon interaction without need of invoking excitonic mechanisms. Our 
work demonstrates that phase-diagrams with competing orders and a superconducting dome are 
also obtained in the framework of the electron-phonon interaction. 
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The interplay between long range charge or magnetic 
order and superconductivity is one of the highly debated 
subject in condensed matter-physics. Doping of an an- 
tiferromagnetic Mott insulator results in high Tc super- 
conductivity in cuprates[l|. In iron pnictides [2] mag- 
netism plays an important role as the parent compound 
is often an antiferromagnetic metal with a large mag- 
netic susceptibility 3] and spin fluctuations are believed 
to survive in the superconducting state 4] . In transition 
metal dichalcogenides such as 2H-NbSe2 or 2II-TaSe2, 
multiband superconductivity [sj and charge density wave 
(CDW) [I coexist. 

Titanium diselenide (lT-TiSe2) belongs to this class 
of materials, as at Tcdw ~ 200K it undergoes a CDW 
instability characterized by a 2 x 2 x 2 real space super- 
structure TiSe2 is not superconducting at low tem- 
perature, but CDW is suppressed and superconductivity 
stabilized either by Cu intercalation (max. Tc — 4.5 K) 
or pressure (max. Tc = 1.8 K) ^. The resulting 
phase diagram looks similar to that of cuprates with the 
difference that the spin-density wave order has been re- 
placed by a CDW oneQ- 

The mechanism at the origin of the CDW and super- 
conducting phases in TiSe2 is currently unknown. The 
similarity with the phase diagram of other correlated 
materials and the semimetallic ^] nature of TiSe2 led 
to the idea [l3| that the transition from the high- 
T semimetallic state to the low-T semimetallic ordered 
stated is of the Overhauser type[llj, namely a Bose- 
Einstein condensation of excitons driven by correlation 
effects. Indeed, as proposed in refs. [l3 . Hsj . in a 
semimetal a charge unbalance is poorly screened due 
to the lack of carriers and the formation of excitons is 
thus favoured. A Jahn- Teller mechanism has also been 
proposed[l3, [lB| to account for CDW in TiSe2 and in 
dichalcogenides in general 116|. This claim is supported 
by thermal diffuse scattering data [l3| showing the occur- 
rence of a soft-phonon mode at zone boundary. However 
the origin of the softening is unclear as no first-principles 



calculations of the TiSe2 phonon dispersion and electron- 
phonon coupling are available. It is then unknown if den- 
sity functional theory (DFT) fails in predicting the CDW 
and superconducting orders. 

In this work we demonstrate that the complete TiSe2 
pressure phase-diagram is reproduced by first principles 
calculations. We show that the disappearing of the CDW 
at finite applied pressure (P) is due to the local chemistry 
around a Ti atom while the pinning of the CDW is de- 
termined by the electron-phonon coupling. Excitonic or 
nesting effects are not the relevant interactions to de- 
scribe the TiSe2 pressure phase-diagram. 

The DFT band structure [1^] and Fermi surfaces of 
TiSe2 at ambient pressure are shown in Fig. ^1 Our 
results are in agreement with previous works |27l . [28| . 
Four bands cross the Fermi level. A narrow Ti-derived 
d-band is almost entirely unoccupied except around the 
L point where it forms elongated and flat hole-pockets (in 
violet). A Se 4p derived band weakly hybridizes with Ti 
states and forms an electron-pocket centered at F having 
the form of a small ellipsoid with larger axis parallel to 
the k^ = plane. Finally the other two bands are formed 
by strongly hybridized Se-p and Ti-d states. Hydrostatic 
pressure weakly affects the band structure except for an 
increase of the Se p band bandwidth. 

The phonon dispersions _18J at pressure P=0,5 and 10 
GPa are shown in Fig. [21 At zero pressure we find that 
the system is dynamically unstable at the L and M points, 
consistent with a 2 x 2 x 1 (M) or a 2 x 2 x 2 (L) real-space 
superstructure in agreement with diffraction data Q . As 
pressure is increased, the CDW disappears at P sa 4 GPa. 
Transport data @ support the occurrence of a supercon- 
ducting state in the region between 2 and 4 GPa. It is un- 
clear if superconductivity coexist or not with CDW order. 
Given that (i) anharmonic effects (neglected here) should 
be important close to the transition to the CDW state 
and should reduce the critical pressure for the occurrence 
of CDW and (ii) the experimental phase diagram at fi- 
nite pressure is only based on transport measurements, 
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FIG. 1: (color online): (Left) TiSe2 electronic structure at P = and 5 GPa and (center and right) P = GPa Fermi surface. 



our results are in excellent agreement with experimental 
finding. Thus DFT correctly describes the occurrence of 
charge density way in TiSe2 at zero and finite pressure. 

To address the microscopic mechanism responsible for 
the occurrence of CDW and its suppression under pres- 
sure we calculate the static bare susceptibility Xo(q) = 



/k+q.i— /k 



and the nesting factor iV/(q) — 



-Wl Sknm '^(ek+q«)'5(ekm), whcrc ekm are the DFT-bands 
measured from the Fermi level. As xo(<l) coincides with 
the electron-phonon contribution to the static phonon 
self-energy assuming constant matrix elements, it repre- 
sents the electron-phonon correction to the phonon fre- 
quencies squared. A minimum in Xo(q) means a soften- 
ing in Wq^. As shown in Fig. [21 Xo(ct) has a minimum 
both at M and L and -/V/(q) peaks at M and L. However 
both quantities are essentially unaffected by pressure be- 
tween and 5 GPa, nesting cannot be responsible for the 
disappearing of the CDW at 4 GPa. 

The mechanism responsible for the disappearance of 
CDW under pressure can be understood by considering 
the force constants C/j(P) at a given pressure P, where 
/, J label the atoms in the 4x4x2 real-space supercell. 
A new set of short-range force constants Cjp{P) is ob- 
tained from Cjj{P) by setting to zero the force constants 
at distances larger than the Ti-Se bond length. Then we 
define: 



Cij{P) - CijiO)-Cf,fiO) + Cff{P) 



(1) 



The force-constants Cij{P) and C/j(0) only differ by 
short-range terms as they both include the same long- 
range behavior calculated at P = 0. From Cij{P) we ob- 
tain new phonon frequencies at any phonon-momentum 
by Fourier interpolation. In Fig. 10 the phonon fre- 
quencies obtained with the aforementioned procedure are 
compared with standard linear-response calculations. As 
it can be seen the two are very similar and, most im- 
portant, the CDW disappears simply by removing the 
short-range force-constants at P=0 and replacing them 
with that at P = 10 GPa. The disappearing of the CDW 
in TiSe2 is not driven by any electronic mechanism or 
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1.57 
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0.84 
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0.62 
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TABLE I: Calculated superconducting properties of TiSe2 us- 
ing McMillan formula (/Lt* =0.1). 



any long-range interaction but it is due to the stiffening 
of the nearest-neighbors Ti-Se force-constants under ap- 
plied pressure. The occurrence of a soft-phonon mode is 
thus inherent to the local environment around the tran- 
sition metal atom in the TiSe2 layer. 

We then consider the superconducting properties of 
TiSe2 at pressures P > 5 GPa. We calculate the electron- 
phonon coupling Xqi, for a phonon mode v with momen- 
tum q and phonon frequency Wq^, namely: 



Wq,7V(0)iVfc 



l5kn,k-Hqml^^(ekn)(5(ek+qm). 



(2) 

El is 5kn,k+qm — 

where V is the Kohn- 



The matrix element in Eq. 
{kn\dV/ i5Mq,y |k -I- qm) / y/2Lu^, 
Sham potential and Uqi, is the amplitude of the phonon 
displacement and iV(0) is the density of states at the 
Fermi level. The critical temperatures obtained with 
McMillan formula and the parameters in table U are in 
qualitative agreement with experiments although in ex- 
periments superconductivity occurs in a narrower range 
of pressures. This discrepancy is probably related to the 
use of a McMillan type of approach to deduce the critical 
temperature and to the neglect of anharmonic effects. 
We calculate the Eliashberg function 

a^F{uj) = -i- Xq^UJqJiuj - UJq^) (3) 



qi/ 



and the integral X{uj) = 2 duj'a^F{uj')/uj' . At P = 5 
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— P=0 GPa but short range force constats at 10 GPa 
-- P=10GPa 
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FIG. 2: (color online) TiSe2 (a) Nesting factor, (b) static 
susceptibility and (c) (d) phonon dispersion as a function hy- 
drostatic pressure. In panel (d) the continuous line phonon 
dispersion is obtained by replacing the nearest neighbours Ti- 
Se force constants in the P = GPa force constant matrix 
with those at P = 10 GPa. 
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FIG. 3: Eliashberg function and integrated Eliashberg func- 
tion A(ix)) at and 5 GPa 



GPa the Eliashberg function is dominated by the soft 
mode [1^ at low energy that contributes to most of the 
coupling. As the pressure is raised, the low energy mode 
hardens (see Fig. [5]) and the average electron-phonon 
coupling is reduced, as shown in Fig. [31 

In the present work we have shown that the finite pres- 
sure phase diagram of TiSe2 can be entirely explained 
by the local chemistry around a Ti atom and by the 
electron-phonon interaction. In more details, the disap- 
pearing of the CDW at finite pressure is determined by 
the stiffening of the Ti-Se nearest- neighbours short-range 
force-constants while the static bare electronic suscepti- 
bility determines the ordering vector at which the CDW 
locks-in. The electron-phonon interaction alone does not 
induce a CDW, as it is clear from the 5 GPa calculation. 

In previous works on other transition metal dichalco- 
genides it was shown that nesting is indeed irrelevant in 
determining CDW order 0-1111 and that DFT is able 
to reproduce the occurrence of CDW in 2H-NbSe2 [111 
and lT-TaS2 In all these works the origin of the 

CDW was attributed to the electron-phonon interaction. 
In our work we show that, although electronic suscepti- 
bility is relevant in determining the ordering vector, the 
local chemistry is the real responsible for the occurrence 
of CDW. More exotic mechanisms such as an Overhauser 
transition Q or excitonic phases ;10'j are unlikely. 

Our work reproduces also the main features of the su- 
perconducting state of TiSe2 as it explains the dome-like 
shape of Tc versus pressure found in experiments Q . Su- 
perconductivity is driven by the softening of the phonon 
mode related to the charge density wave. As the mode 
is very soft, substantial anharmonic effects are to be ex- 
pected. The inclusion of these effects would harden the 
phonon frequency of the soft mode and shift the dome re- 
gion in which superconductivity takes place to lower pres- 
sure in even better agreement with experimental data. 
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Several authors f?!, ^ pointed out the strong similarity 
of the phase diagram of TiSe2 with cuprates and sug- 
gested a non-conventional pairing mechanism for TiSe2. 
Our work shows that this similarity is only apparent 
as here superconductivity is phonon-mediated. Further- 
more our work demonstrates that exotic phase diagrams 
including dome-like superconducting regions, as those oc- 



curring in Cuprates, Pnictides or Cobaltates can also 
arise from an electron-phonon mechanism. Thus, the ap- 
parent similarity of these phase diagrams to the cuprate 
case cannot be considered a fingerprint of electronic cor- 
relation. 

Calculations were performed at the IDRIS supercom- 
puting center (project 081202). 
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